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1 Introduction 



Many efforts in contemporary physics are devoted to the study of relativistic 
quantum systems in non-equilibrium states, like e.g. the expansion of hot, com- 
pressed nuclear matter. Unfortunately, the combination of relativistic quantum 
field theory and statistical mechanics required for the theoretical description of 
such systems, i.e., relativistic transport theory, is far from well understood and 
still subject to discussion. As a result, most publications in this field either start 
from field theory and end up with some transport equation that can hardly be 
solved - or solve a transport equation, which has lost much of its touch with 
quantum field theory. 

The present paper attempts to bridge this gap through a different approach: 
The response of a quantum system to gradients in temperature and chemical 
potential is calculated directly, i.e., without the detour of transport theory. To 
this end a formalism is introduced, which incorporates the space-time depen- 
dence of temperature and chemical potential into the quantization rules for the 
interacting quantum fields. 

The paper is organized as follows: A brief introduction to the formalism of 
Thermo Field Dynamics is followed by its application to a free charged scalar 
field in global as well as local equilibrium states. The ideas are then extended 
to the interacting field, and used to calculate the full two-point function for 
spatially inhomogeneous states approximately. 

In section 5, the formal results are then applied to calculate the currents of 
the system, which arise due to the spatial inhomogeneity. Kinetic (or transport) 
coefficients are obtained in terms of the spectral function of the interacting field, 
and compared to formal results obtained with different methods. Numerical 
results are presented for a simple approximation to pions interacting with a 
nuclear medium. 

An important aspect of finite temperature field theory is the existence of 
two different, mutually commuting representations of the canonical commuta- 
tion relations: The Hilbert space of the quantum description is "doubled" with 
respect to the zero temperature case Q. For a description in terms of Green's 
functions, this requires to use 2x2 matrix valued propagators. The two best de- 
veloped realizations of such a formalism are Thermo Field Dynamics (TFD)|Q 
and the Closed-Time-Path method || . 

While the vacuum state of quantum field theory is invariant under Lorentz 
transformations, matter states have, in general, less symmetry. This leads to a 
modification of the particle spectrum needed for the description of such a state 
[|J. Hence, a consistent quantum description requires elementary excitations 
with a continuous mass spectrum, rather than physical quasi-particles of infinite 
lifetime §. 

Formally this is expressed as spectral function of "particles" which deviates 
from a (5-function, i.e., which is more than a mass-shell constraint. In several 
papers it was shown recently, that in the framework of Thermo Field Dynamics a 
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physical interpretation of the 2x2 matrix structure is obtained also for systems 
with such a nontrivial single-particle spectrum J|, [| fjj. In this sense, TFD is 
superior to the Closed Time-Path method. 

In "ordinary" quantum mechanics, a thermal state of a quantum system is 
described by a statistical operator (or density matrix) W, and the measurement 
of an observable will yield the average 

v ' 7 Ti[W] w 

where the trace is taken over the Hilbert space of the quantum system and 
£ is the hermitean operator associated with the observable. In Thermo Field 
Dynamics (TFD), the calculation of this trace is simplified to the calculation of 
a matrix element 

\ _ jW L \\£(t 7 x) \\W*} 
^' j / {W L \\W*)) ' [) 

with "left" and "right" statistical state defined in terms of the two different 
commuting representations (see refs. ^ ^| for details). 

If gradients of the temperature etc. are present in the state, the hamiltonian 
of the quantum system, which generates its time evolution, contains a part 
which is due to the spatial inhomogeneity. The above statistical average in the 
course of time evolution will be modified by the gradients present in the system. 

If the space-dependence of the temperature and chemical potential is suffi- 
ciently weak, the response of an observable £ to a gradient is a linear function 
of the perturbation. We thus express this response 5 (£(t, x)) up to first order 
in Hi (the gradient part of the hamiltonian) as 



5(£(t,x)j =iJdT{j£(t-T,x),Hi\ J . (3) 

to 

If the commutator can be calculated, one can use it to obtain the transport 
coefficients of the system, i.e. the coefficients relating the response to the per- 
turbation. The crucial point however is the commutator: Its expectation value 
is generally calculated as 



£{t-r,x),H^j = Jd\(^- t £{t-T-x\,x)H^j , (4) 
o 

i.e., as the time derivative of a correlation function with complex time arguments 
[|[ [)| . The above expression is usually referred to as the Kubo formula [[ToJ . 

The analytical structure of a correlation function in the complex time plane 
is related to the properties of its Fourier transform in the complex energy plane. 
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However, the latter is plagued by pathologies of perturbation theory which are 
especially severe when considering nuclear phenomena |ll], [l2| . The most natu- 
ral way to overcome this problem is to use TFD in a direct computation of the 
commutator expectation value of eqn. (0). 

We therefore consider a complex, scalar boson held describing spinless charged 
excitations in a statistical system not too far from equilibrium. In the spirit of 
the first remark, one could think of this field as describing positive and negative 
pions in nuclear matter. 

According to the reasoning above, the thermal boson field is then described 
by two field operators <j) x , 4> x and their adjoints </>*, 4>*, with canonical commu- 
tation relations 



4>{t,x),d t 4t*{t,x') = 



= iS 3 (x-x') 
- i8 3 (x-x') 



(5) 



but commuting with each other. They are combined in a statistical doublet 



and the 2x2 matrix valued propagator of this doublet is 



(6) 



(a 



= -K(^ L ||T 



\\W R ))/{W L \\W R )) 



iW L \\T [^ X ,]\\W R )) {W L \\T 
iW L \[T fe^l 1^1 iW L \\T 
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\W R \ 

\W R )) I lw L \\w R )) 



7) 



The free as well as the interacting scalar field can be expand into momentum 
eigenmodes 



d 3 k 



t u\ n -ikx , / ,\ Akx 
Ofc_(<)e +a k +{t)e 



-^=(4_(t)e ifc -+S fc+W 



ikx 



(8) 



h is the three- momentum of the modes, therefore in this notation at (t) creates 
a negatively charged excitation with momentum k, while a&+(t) annihilates a 
positive charge. Henceforth the two different charges are distinguished by an 
additional index / = ± whenever possible. 

For the free case the commutation relations of the a-operators at different 
times are simple, while they are unknown for the interacting fields. Note, that 
in the above expansion we have used a different energy normalization factor 
than usual. While this might be a little weird for free fields, it simplifies the 
calculations for the interacting case. 
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2 Free charged boson field 



The free complex scalar quantum fields are expand as specified in eqn. (j^) . The 
operators have commutation relations 



a k i(t),a\ lv (t) 
aki{t),a} k , v {t) 



1 



1 



S w S 3 (k - k') 
6 W S 3 (k - k') , 



(9) 



for I = ±. All other commutators vanish, and w is the free on-shell energy 
depending on momentum and mass m of the physical particles as 



L0 k 



k 2 



(10) 



The free physical states created and annihilated by the above operators are 
stable. However, they cannot be used as input for perturbative calculations, 
since an infinitely small interaction gives them a finite lifetime ]lj . 

The stable (albeit unobservable) quasi-particle states are related to the phys- 
ical states through a Bogoliubov transformation. The hamiltonian of the system 
is invariant under this transformation, if the system is in global equilibrium. 
Hence, in this case the system has a global Bogoliubov symmetry. The trans- 
formation can be written in matrix form as 



««(*) 
-aid(t) 



-c 



kl 



€kl 



B(m(k))- 



l 



(11) 



The ^-operators create and annihilate stable momentum eigenmodes, and the 
symbol ^ has been chosen to indicate, that this operator is not the hermitean 
adjoint of £. Rather, these operators are thermal conjugate to each other, i.e., 
they annihilate the thermal ground states according to 



Zm\\w r )) = q , li\\w R )) = o , lw L \\& = Q , ((w L ||££ = o vfc,; = ± 



and have commutation relations 



= 6 u ,6 3 (k-k') 
= 5 w 5 3 (k-k') 



(12) 



(13) 



(all other commutators vanish). Using these rules, any matrix element of field 
operators can be expressed through the elements of the Bogoliubov matrix. 
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The Bogoliubov matrix can be parameterized in various ways, leading to 
unitary inequivalent representations of the thermal quasi-particles ]2|, Q |i"3[ . 
The most useful parameterization however is obtained as linear function of a 
single parameter n that resembles the phase-space occupation factor, 

S(n) = (^ • (14) 

In general, the n-parameter can be different for any mode of the system. The 
modes of the free scalar field are characterized by their momentum, hence n = 
n{k). The field is complex, implying that the bosons carry a charge and thus 
are not their own anti-particles. Hence, for the two species of bosons, the n- 
parameters can be different, n + ^ n_. If the charge they carry is a conserved 
quantity, and the field is in thermal equilibrium, we can chose Bose-Einstein 
functions as 

n± (Jfe) = /±(fc) ? - = __L /±(fe) = e-/»<«*=w) . (15) 

v ' l-/ ± (fe) e^Kw)-l v ' v ' 

Note, that in principle a generalized equilibrium state with (3 = (3(k) could be 
used as well. 

It is an easy task to calculate the propagator from the above considerations, 
we refer to the existing literature or the next section for its explicit form. More 
instructive for our purpose is the case, where one uses a Bogoliubov transfor- 
mation that mixes the momenta as [fTil 



—\uj a t 



T 



a ^}A = -tL= fd 3 q[ ^ l I B^k)^ . (16) 



This momentum mixing occurs, when one makes the Bogoliubov symmetry 
mentioned before a local symmetry. Such a process is called the gauging of 
the symmetry, and in a quite natural way introduces a couping to external 
parameters into the model |fl5| . The generalized Bogoliubov matrix is expressed 
as 

It obeys the relations 

(B^fok))* Bi{q,k') = 5 3 (k-q)5 3 (q-k')dmg(l,l) 

+ T r 3 N l (k,k')(s 3 (k-q)-8 3 (q-k') 
fd 3 q (&r\q,k)y Bi(q,k') = 5 3 (k - k') diag(l, 1) , (18) 
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with the Pauli matrix t% — diag(l, — 1) and 

n=( \ J) . (19) 

The occupation number density parameter appearing here is the Fourier trans- 
form of a space-local quantity 

JV,(9 '* !)= (2^ f d3ze ~ kq ~ k}Z M(<l + k)/2,z) , (20) 

and we specify that the n; (k, z) are real functions. This implies, that N*(q, k) = 
Ni(k, q). Moreover, the generalized Bogoliubov matrix from eqn. ( |l7j ) is related 
to the one from (|l4|) as 



Bi{q, k) = JSze-W- k ) z B(m((q + k)/2, z)) . 



(21) 



As one can check within a few lines, the above transformation preserves the 
canonical commutation relations. However, the time evolution is complicated 
by the simple fact that the phase factors in eqn. ( |l6|) carry the momentum index 
of the ^-operators. Indeed, one obtains for the time derivative of the a-operators 

d ( a kl (t) \ / a k i(t) 



+ [d 3 q , E N t (k, q) K - Uq ) T r 3 ( ) .(22) 

A similar relation is obtained for the adjoint operators. Note, that the additional 
term in this time derivative vanishes, if N(q, k) is proportional to 6 3 (q — k). 
This is the case, when n(k, Z) does not depend on the space coordinate Z, i.e., 
when the system is translationally invariant. 

We can now ask, what kind of Hamilton operator H would give a time 
derivative according to eqn. (E2J) if used in the Heisenberg equation 



.d_ ( a H (t) 
l dt \ al(t) 



au (t) 



(23) 



#o = E/^(_SjJ { S ) . (-) 



Clearly, such a hamiltonian has two parts: one "bare" giving the unperturbed 
time evolution 

and a part that vanishes for homogeneous systems: 

a J$)'t*,*,m ($> ) , t») 
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with 

Hi(k,q) =Ni(k,q)2^uT q (w* - w g ) . (26) 

Note, that the dimension of the numerical factors in ( pi[ ) and ( |25| ) is mass 2 , due 
to the energy normalization factors contained in the a-operators, cf. eqn. (^|). 
The further discussion of this modified hamiltonian we postpone until the above 
expression has been generalized to interacting systems. 



3 Interacting fields in global equilibrium 

The procedure to derive thermal Green's functions for interacting fields has been 
outlined in refs. H 0], and thus will not be repeated here in full detail. It is 
based on the expansion of the interacting field into modes with definite energy 
and momentum according to ref . |l| . Their superposition with a positive weight 
function then gives the fully interacting field. 

Such an expansion is necessary, because at finite density and temperature the 
irreducible representations of the space-time symmetry group are characterized 
by two continuous parameters Since these fields contain the full interaction 
of the system, the commutation relations of the fields with definite energy and 
momentum are not known in general. However, for the task of calculating 
two-point functions or bilinear expectation values, one needs to know only the 
expectation value of such commutators - and this information can be absorbed 
into the weight function. 

In other words, for the calculation of bilinear expectation values of interact- 
ing fields it is sufficient to consider the modes of definite energy and momentum 
as generalized free fields [Q . The full information about the single-particle spec- 
trum of the theory is contained in the weight functions. 

The final step in deriving full propagators at finite temperature then consists 
in applying the thermal quasi-particle picture to each of these modes. We thus 
write for the operators of the field expansion 

~HS ) = JdE P T(E,k) (B(n±(E, fc))) -1 ( |f fc± ) e-- 

a k±\ l ) / \ S£fc± / 



= JdEp^ 2 (E,k) I i Ek ± 1 B(n ± (E,k))e^. (27) 

This Bogoliubov transformation for interacting systems defines stable, albeit 
non-observable quasi-particles || ||, [?J , created and annihilated by the ^-operators. 
The ^-operators have different commutation relations, cf. ref. |Q 

Zeu^m] = S u ,6(E-E')5 3 (k-k') 
hkiM'k'A = S u ,6(E-E')S 3 (k-k') (28) 



tt 



T 
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(all other commutators vanish), but the thermal state conditions remain similar 
to fill) 



£ Ekl \\W R ))=0 , £ Ekl \\W R }=0 , ((W L \\^ kl =0 , lW L \\t$ kl =0 VE,k,l = ±. 

(29) 

With these rules, all bilinear expectation values can be calculated. Matrix 
elements of higher powers of interacting field operators however are obtained 
only approximately, they have a perturbative expansion in terms of the weight 
functions p t (E, k) [|. 

n is a Bose-Einstein distribution function, which for special case of a global 
equilibrium depends only on the energy of the modes: 



, f±(E) 



(30) 



In our notation, the weight functions pi(E,k) have support only for positive 
energies, and their integrals are 



oo 

JdEE Pl (E,k) = I 

o 

oo 

JdE Pl (E,k) = Z kl . 



(31) 



Note, that the existence of this spectral decomposition is only guaranteed in 
case the system is space-time translation invariant, i.e., if it is in a thermal 
equilibrium state. It follows from these definitions, that 



aki (t),a\, v (t) 



Z kl 5 W 6 3 {k-k') 



,j(t),ojL(*)l = ZkiS u ,s 3 (k-k') 



(32) 



in generalization of eqn. (||). Collecting these definitions, the full bosonic two- 
point Green's function, or propagator, at finite temperature is derived as || M 



D^(k ,k) = 
dEp+(E,k) (Bin+iE^k)))- 1 

dEp_{E,k) T 3 B T {n_(E,k)) 



1 



fc - E + ie 



1 



ko — E — ie 



V 



1 

E + ie 



B(n + (E, k))r 3 
(B T (n^(E,k)))m 
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In the above equation, the terms propagating particle and anti-particle states 



have been kept separately. Eqn. (30) can be continued to negative energy 
arguments, and 

n-(E,k) = -(l + n+(-E,-k)). (34) 
Setting n+(E, k) = n(E) in the following, the above propagator simplifies to 

L» (Qf>) (fc ,fc) = J dE A(E, k) x 

— OO 

1 



(Bl/zlD))- 1 | Ko-E + ie x | B[ll{E]]Ti , :t] 

fco — E — ie 

A(E, k) is the spectral function of the boson field, 

A(E, k) = p+(E, k)Q(E) - p-(-E, -k)Q(-E) , (36) 
and the limit of free particles with mass m is recovered when 

A(E, k) — ► sign(£J) S(E 2 - k 2 - m 2 ) = sign(E) 5{E 2 - Lu 2 k ) . (37) 

4 Interacting fields with local Bogoliubov sym- 
metry 

As was stated above, the existence of a spectral decomposition is bound to the 
equilibrium property of the system, i.e., to its space-time translation invariance 
[||. We can expect however, that close to equilibrium the field properties do 
not change very much. Thus, the considerations from the second section of the 
present paper can be carried over to the interacting case: We want to study 
the quantum system under the influence of small gradients in the distribution 
function n, but with the approximation of a space-time independent spectral 
weight function pi(E, k). In generalization of eqn. (ttq): 



= SdEjd*qpy\E,(q+k)/2) (Bf\E,q,k)) { ^ j e— 

(_~«(t))' = ] dE Sd 3 <lP 1 AE,(q + k)/2) ( ^ S ; (S,qf,fe)e^{38) 

The transformation matrix here has the same functional form as in eqn. (|l7|), 
but now its occupation number density parameter also depends on the energy 
variable 

Nl (E,q,k) = ^ J^Ze-^-^ z m(E,(q + k)/2,z). (39) 
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The propagator matrix according to (0) is calculated as 

oc . 

D^(t,X;t',x')=-i JdE j * k f 2 «f k P?(E, (k + q)/2) fif(E, (q + fe')/2) 
o 



-e(f - 1) 

i IdE 



d3k ^f k ' P^(E, (k + q)/2) p^(E, (q + k')/2) 



r 3 B T (E,q,k) ^ ~^ _ Q{t _ tl) ) (gE(E,q,k'j) ~ ^-^-^+^(40) 

Similar to the equilibrium case and to the time-dependent case one can perform 
the energy integration because of the linearity of this expression in n. This allows 
to write the full propagator in terms of Bogoliubov matrices, which depend on 
space-local distribution functions. The principles of this method are outlined in 
ref. ||, and will not be discussed in the present paper. 

As the last point in considering the above full Green function, consider the 
imaginary part of the retarded propagator obtained as 

lm(D R (t,X;t',x')) = lm(^D^(t,X;t',x')- D^(t,X;t',x') 

rdEd 3 k 



The retarded and advanced propagator hence only depend on the coordinate 
differences, and their imaginary part still gives the spectral function from eqn. 
(j36l). This is the only approximation made so far, but it can be shown easily by 
solving the matrix-valued Schwinger-Dyson equation in coordinate representa- 
tion, that corrections to the above spectral functions are of second order in the 
gradients. 

In generalization of eqn. (|22), one thus obtains for the time derivative 
d ( a a (t) \ ( a k i(t) 



X 



E P ^(E,k)-E P y 2 (E,q))T r 3 ( ) .,42, 
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Here, the energy ft is defined as 

oo 

Sl u = ^- fdEE Pl (E,k) = -j—, (43) 



and Zki as in eqn. ([H]). Obviously, as for the free field, the additional term in 
this time evolution vanishes for translationally invariant n, i.e., when N(E, k, q) 
is proportional to S 3 (k—q). It has a generator Hi, which has the same functional 
form as in eqn. (|25|), but the kernel reads 

TUM = JdE P -I {E z l k z + q q)/2) m(E,k,q) ( Ep y\E,k) -E p y\E,q] 

(44) 

The perturbation of the system has a very special form given in eqn. (|25|). If it 
is expressed in terms of the ^-operators, one may use 

B(E,q,4)(\ j)(B- 1 (S,fc,fc'))*=(° ~J ) S 3 (q-q')S 3 (k-k') 

(45) 

to see that only the combination £tt£ft survives: 
00 

Hi = JdEdE' Jd 3 kd 3 k' p 1 / 2 (E,k)p 1 / 2 (E',k') Ui(k,k') ZtkM'vi ■ 



i=± 



(46) 

This simple expression for the gradient part of the hamiltonian can be inserted 
into eqn. (|J). Note, that the total hamiltonian is diagonal in the ^-operators, 
but not if expressed in terms of the physical particle operators a, a*. 



5 Transport properties of an interacting scalar 
field 

Observables are expressed as functionals of the interacting fields - and therefore 
as functionals of the operators a, (V. Hence, although the state we consider is 
stationary in terms of the basis defined by the ^-operators, momentum mixing 
as introduced above will result in a non-trivial time evolution of physical quan- 
tities, according to the Heisenberg equation (p3|). Thus, equation (Q) can be 
used to calculate relevant contributions to system parameters due to the spatial 
dependence of n(E, k, x). 

Two quantities we are interested in are the conserved particle current oper- 
ator of the complex scalar field and its expectation value, 

5(t,aj)=i(^V0 x -0 x V^) j(t,x)=s(j(t,x)) (47) 
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and the conserved energy current 

%,cc) = (^*V0 !C + 9^ x V</»*) E(t,x) = 6(E(t,x)j . (48) 

Inserting these into eqn. (|^) then yields an expression for the commutator, 
which contains four field operators instead of two. The expectation value of this 
commutator therefore cannot be expressed completely by the interacting two- 
point functions we have calculated. For these, it had been explicitly assumed 
above, that the fields can be written as generalized free fields. 

However, the four-field expectation value has a diagrammatic expansion in 
terms of the propagators we derived - and the lowest order term of this expansion 
already contains the full single-particle spectrum. We expect, that this lowest 
order contributes the dominant part of the system's response. 

Using the commutation rules for the generalized free fields, e.g. the Bogoli- 
ubov transformation and (p8f), therefore amounts to the neglection of higher 
correlations. These however can be incorporated through the appropriate Feyn- 
man rules [[!} [ill . 

If we insert the expansion of the fields according to (||), all kinds of binary 
products of the et-operators appear in these operator expressions. However, in 
the present formulation, out of 4 x 8 combinations of a-operators contained in 
a naive calculation of the commutator (||), only two survive due to the special 
form © of Hi. These are 

ahit - T ) a k'i(t - T),^ E , q , l £} Eql 

= pl /2 (E, q) P ] /2 (E', q 1 ) S 3 (k - q) 8 3 {k! - q 1 ) ^B-E>){t-r) 

9i4;(t-r) a k >i(t-T),^ E , qll l 

= iEp] /2 (E, q) p¥ 2 (E', q') 5 3 (k - q) S 3 (k' - q') e i( - E - E '^- T \49) 

for I — ±, i.e., both charge species. 

We can then assemble the results for the currents, by performing the time 
integration as if t — > oo and to — * — oo. Alternatively, one could first process 
the momentum integrations and then pick out the real part of the current. One 
obtains 

i l (t,aj) = -27ri JdE J'^^p l {E,k)p l (E,k l )e^ k - k '^ (k + k') H^k') 
o 

(50) 
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for each species. The energy current is formally quite similar, 

CO f 

E l (t,x) = -2m JdE 1^1^- pi(E,k) Pl (E,k')e^ k - k '^ E^k + k') H t (k,k') 
o 

(51) 

and the total currents are for both cases 

j(t,x) = j + (t,x)-j_(t,x) 

E(t,x) = E + (t,x) + E_(t,x) . (52) 

In principle the above expressions can be calculated, when the spectral functions 
and the space-dependence of n are given. However, it is more instructive to 
perform a gradient expansion for the occupation number density n. For this 
we use the abbreviation Q = (k + k )/2, and expand all momenta around this 
value, e.g. 

p 1 l /2 (E,Q)(Ep 1 l /2 (E,k)-Ep 1 l /2 (E,k / )) « (k-k')^^L Pl (E,Q). (53) 

The factor (k — k') can be taken out of the integrations, as —id/dx = — iV x - 
Other contributions are of zeroth and second order in this derivative, e.g. 

Pl (E, k) Pl (E, k') = pf(E, Q)+0 ((fc - k'f) . (54) 

Therefore, the i-th vector component of the l-charged currents generated by the 
inhomogencity of the system is 



jf{t,x) = 27 r/ggg^ JdE( Pl (E,Q)) 



2 

X 



dE , E , fdni(E',Q,x)dp l (E',Q) 



Ef\t,x) =2nf0lQ^- JdEE (pi(E, Q)) 



, 2 

X 



J \ dxO) d Qh) j 

Before this is evaluated further, note that the expression in the curly brackets 
is nothing but the first term of a full gradient expansion of the product of n and 
p, i.e., the Poisson bracket. It is therefore clear, that the above expression is 
useful also for situations, where the space-dependence of p cannot be neglected: 
the Poisson bracket then also has a contribution with a space-derivative acting 
on p and a momentum derivative acting on n. The above expression is also 
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interesting in view of the original equation (^) , since the Poisson bracket is the 
analogon of the commutator, i.e., it contributes a factor U to the the current. 

To proceed, some remarks on the derivatives occuring in the Poisson bracket 
are necessary. It was assumed, that the system is in a local equilibrium state, 
close to global equilibrium such that the gradients are small and currents are 
only due to these gradients. This implies, that the spectral function depends 
only on Q = \Q\, and thus the momentum derivative gives 

U) 



dpi(E',Q) = Q U) dpi(E',Q) 
qqU) Q dQ 



(56) 

The currents from eqn. (^5|) thus are symmetric in the indices i and j. 

Using as local equilibrium distribution function instead of eqn. ( |30| ) the 
expression 

the derivative with respect to X has two components, 



dxtt)'"^ , ~ 6, "' J 4sinh 2 (/3(£;T^)/2) \ T T T 

(58) 

This derivative however has to be taken with care: In eqns. (p35|), the X- 
derivative was shifted through the integrals as if they were all properly defined. 
This requires, in the absence of a Bose condensate, that the spectral function 
vanishes linearly at E — fi. A discussion of this requirement will be attempted 
in a separate publication, here we will take it for granted. 

However, this violates the assumption of a spectral function independent of 
the gradients, since the ^-derivative then also acts on p(E,k). To avoid this 
complication, we define the integral over the distribution function as a principal 
value integral, i.e., it is only defined in the sense of first calculating 



I Q l = JdEp(E,Q)n(E,Q,x) = lim J + J 

L u+e 



A*— « oo 



dEp(E,Q)n(E,Q,x) 



(59) 

and then taking the CC-derivative. The same decomposition into T and p- 
gradient terms then holds for the currents, i.e., with transport coefficients Lij 
we obtain |j] 

VT 

3i(t,x) =-L n Vn -L 12 — 

VT 

Ei(t,x) =-L 21 Vp -L 22 — . (60) 
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Note, that in our picture L12 and L21 can be different, i.e. the Onsager relation 
Lij = Lji is not necessarily fulfilled. The reason for this is clear: in contrast 
to ordinary many-body quantum physics, the present formulation exhibits dis- 
sipation already on the tree-graph level jlj] . In other words, the formulation of 
an interacting theory with continuous spectral functions is not micro-reversible, 
physical states have a finite lifetime. 

In summary of the above reasoning, we can thus write the kinetic coefficients 
with the abbreviations 

X Q i = JdE (pi(E,Q)f 

Y Q i = JdE E ( Pl (E,Q)) 2 , (61) 

plus the definitions of Iqi from eqn. (|59|), Zqi from eqn. (31) and Q^Q^ — > 
Q 2 /3 in very compact form as 

T fd 3 Q Q ^ . m X Q i d d 

L*Y\ = —2n / — > sign / — — — — — Ini 

J 2tt 3 6 ^ ' ' Z 2 n , dfi dQ Q 



(27r )3 6 ^ w Z 2 Ql dn dQ 
Ll2 = ^7(2^6 I^^Z^JdfdQ 1 ' 



fd 3 Q Q^X^Jlj^ 



~2n 



d 3 Q Q^Y Ql d d 
(2^)3 6 ^ZlJdfdQ 1 ^- 



By subtracting the convective part from the energy current, the thermal con- 
ductivity for the interacting scalar field is obtained as 

L11 = d is the diffusion coefficient, and in analogy to electromagnetic plasmas 
we call the quantity 

a =TL^ (64) 

the thermo-force on the charged bosons. 

The above formal results can be compared to different calculations on a 
numerical as well as a formal level. To begin with the latter, we insert a simple 
spectral function for the boson field at zero chemical potential, 



2 E lk 

tt {E^-niy + AE^i 



p+ (E, k) = p_ (E, k) = — ^ — 2 ; j^- 2 (65) 
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with 0% = Cf. + 7j! and jk -C £fc- Such a spectral function can be considered 
the lowest order approximation to a system with nontrivial self energy function 

In ref. JL8J , this prescription was applied to the pion dispersion relation in 
nuclear matter by rewriting the inverse retarded propagator for the pions with 
real energy E as 



E 2 - El(k) - U R (E, k) = {E- {e k - i 7fe )) (E + (e k + hk)) 



(66) 



With this on-shell approximation the integrals are 

1 



dEpi(E,k) 



2ire k 
1 

2iT 



7T 

2 ~ 

Ik 



arctan 



ll 



and 



X k i 
Yja 



1 



7T 2 £k 



^4 



•0( 7 ) 



-o( 7 4 ; 



(67) 



^+0(7) 



(68) 



In the absence of a chemical potential, L12 and L21 are zero, but diffusion 
coefficient and thermal conductivity are to lowest order in the width 



d = -± 
T 



k 1 d / w 



A 



1 

T 



(2tt) 3 3 7fe 9fc 

(/ fc fcefc d ^ + n ( £fe )) 



(2tt) 3 3 7fe 5A; 



(69) 



Here, n(efc) is the local Bose function taken at energy e k - Apart from the 
momentum factors, a similar representation for A was obtained in ref. The 
difference can be attributed to the fact that in this reference a hydrodynamical 
picture was assumed together with f3jk -C 1. We thus make a comparision in 
the high temperature limit, where one obtains 



lim A = — T 

T^oo 7 



(70) 



Here, 7 is some momentum averaged width and A a numerical factor. The same 
limiting behaviour is obtained in |||. Also other calculations in quantum field 
theory give qualitatively similar results: The bulk viscosity of matter should 
rise with the same temperature dependence as the thermal conductivity, and 
in [[l9] (and references quoted there) was found to rise as T 3 in a quark-gluon 
plasma. 
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The above expressions are also in accordance with naive expectations: The 
thermal conductivity and diffusion coefficient diverge, if the width of the parti- 
cles is reduced, i.e., if the interaction is removed. In this case the time needed 
for the relaxation of a temperature disturbance is infinite. 

6 Application to pion propagation in nuclear mat- 
ter 

In the next step one has to look at numerical results. In the spirit of the 
introductory remarks, pions were chosen for this example. However, to obtain 
meaningful results one has to go beyond the simple approximation discussed 
above. Instead of the on-shell approximation of eqn. ([36]), the so-called A- hole 
model at finite temperature is used |Ts| ]. 

In this model, pions couple to excitations in a gas of (non-relativistic) nu- 
cleons and A 33 resonances of width V. Although a numerical calculation of 
the pionic spectral function in this system is not very difficult, for the present 
calculation we use an analytical approximation 

, k 2 C TEuo A 



n {E 2 - w\ ) 2 [E 2 - uj 2 _) 2 + T 2 E 2 {E 2 - E 2 ) 2 ' 

(71) 



E 2 = yk 2 + rn 2 is the free on-shell energy of the pion, uj a and lu± are momen- 
tum dependent, 



uj a = E A (k) - M N = J k 2 + Ml - M N 



En a = V (wa + g' C) 



I = \ [E NA + (r/2) 2 + El ± J (E 2 NA + (T/2) 2 Elf + 4fc 2 Co; A ^) 



and C is the temperature dependent effective coupling constant (see table) 



£ 7T 

JNA 


9' 


m w Mn 


M A 


r 


2 


0.5 


0.14 GeV 0.938 GcV 


1.232 GcV 


0.12 GeV 



Table 1: Coupling constants and masses used in the calculations of this work. 
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with partial densities p° N and °f nucleons and A's. The partial densities at 
a given total baryon density ps = p° N + pA are obtained by a self-consistent 
calculation of the chemical potential fiN = P-A of the fermionic sector. 

The above spectral function is obtained with two approximations: First, its 
is correct only to lowest order in the baryon density (quasistatic approximation) , 
and second it is correct only to lowest order in r/M^ (asymptotic expansion of 
the A-spectral function for constant T). These assumptions yield a polarization 
tensor of the A-nucleon gas 

n "< E - fe > = fe2c WTW^T ,74) 

which is then subject to a Migdal (or g') correction (see refs. [|l8| |2(], ^l] for 
more details on this model) . Figure [j] is a logarithmic contour plot of the above 
spectral function at a total baryon density of pb=1.69 nuclear matter density 
and temperature T = O.lGeV. It exhibits two branches, corresponding to the 
mixing of a coherent A-hole excitation and the pion. At low momenta, the 
strength of the excitations lies predominantly on the lower branch, while at 
higher momenta the pionic strength has changed to the upper branch. 

In figure @, the product of thermal conductivity A and temperature T which 
one obtains with this spectral function is plotted as function of temperature at 
different densities. Figure |^ is a similar representation of the diffusion coefficient 
d. 

A comparison of these values to experimental data is beyond reach for the 
time being, we thus have to restrict the comparison to other calculations. How- 
ever, only very few calculations based on field theoretical methods are available. 
We have already mentioned the agreement with ref. Q on a formal level, if a 
simple approximation for the spectral function is inserted (eqn. (pa)). Although 
the structure of the more sophisticated spectral function from ( |7l[ ) is quite dif- 
ferent, the high temperature behaviour at the example density of pi, = 1.69x 
nuclear matter density inferred from figure ^| is 

A w 9897 GeV/fm 2 x [T/GeV] 3 ' 70 , (75) 

i.e., it rises even faster than estimated by eqn. wW. This is due to the decreasing 
width of the pion with temperature in the A-hole model jl8| - a fact, which 
was not included in eqn. (70). At lower temperatures, the thermal conductivity 
rises slower with temperature than the asymptotic expression. 

The density dependence of the coefficients A and d is quite small, with the 
general tendency to have a "stiffer" temperature dependence at lower baryon 
density. This is again consistent with the fact, that the pion in the A-hole model 
becomes a free particle without the baryonic background. 

The exploratory calculation of the present paper was done with a constant 
width r for the A-resonance. Even the vacuum spectral function of the A is 
not completely described by such an ansatz, much less its medium dependence 
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[ pl| p2[ . The width in matter at finite temperature (obtainable by self-consistent 
calculations currently in progress) might be a factor 2-3 times higher, hence 
according to eqn. (^9|) the thermal conductivity should be lower due to the 
medium dependence of the A-width, and its temperature dependence flatter. 

For completeness, also variational and hydrodynamical calculations of the 
pionic thermal conductivity in the baryon-free regime have to be mentioned 
[ p3[ Q (also plotted in figure ||). In these, the thermal conductivity for T ~ 
0.1 GeV is much lower than obtained in the present paper: Beyond T « 60 
MeV it rises only quadratically with temperature (see thin lines in figure || for 
the corresponding slopes). The results have been obtained by using vacuum 
pion-pion scattering data as input, and thus the interaction in these models is 
completely different from the A-hole model. Hence, a direct comparision of the 
results is not justified. 

Nevertheless, a short remark on the difference is necessary. Within the 
hydrodynamic picture a medium effect, and thus also the Independence of the 
thermal conductivity, can be explained by the rise of the effective pion-pion 
scattering cross-section with temperature p5[ ]. It was already stated, that the 
medium broadening of the pion might be underestimated in the present paper. 
However, even when reducing the present results by a factor 2-3 at higher 
temperatures, we still find that the thermal conductivity of the interacting pions 
in presence of the hot baryonic background is higher than in refs. |23, GM. 



7 Summary and Conclusions 

In the present paper, Thermo Field Dynamics (TFD) for inhomogeneous sys- 
tems |^|, [lJJ is generalized to quantum fields with a continuous single-particle 
mass spectrum jjj. Following a procedure outlined in refs. [^], Q] , the spectral 
and statistical content of the full two-point function are separated in a well de- 
fined approximation. This approximation, i.e., assuming that a spectral function 
still is a meaningful concept, limits the results to situations not too far from 
equilibrium. However, corrections to this approximation are of second order in 
the gradients. Furthermore a general formulation in coordinate space is easily 
achieved (5[ |(| , which also takes temporal inhomogeneity into account and thus 
does not trade simplicity for covariance. 

By requiring the thermal Bogoliubov symmetry to hold locally in space, a 
modification of the systems' hamiltonian was obtained. It is well known, that 
this gauging of a global symmetry introduces a minimal coupling to an external 
classical field, here: A "temperature field" and a "chemical potential field" p!q |. 

In a system close to equilibrium, this coupling was used to calculate trans- 
port (or kinetic) coefficients of the interacting quantum field. In contrast to 
previous derivations neither imaginary time arguments nor perturbation 
theory in powers of a coupling constant were used for this purpose. Expressions 
for thermal conductivity and the diffusion coefficient were given in terms of the 
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equilibrium spectral function of the quantum field. They are dominated by two 
factors: A Poisson bracket, which stems from a gradient expansion (cf. eqn. 
(p5[)), and proportionality to the inverse width of the spectral distribution, cf. 
eqn. @. 

The Poisson bracket contributes a factor h, while in a naive counting the 
inverse width contributes a factor 1/h to the transport coefficient. This can- 
cellation seems to be the main source for the difficulties one has in obtaining 
kinetic coefficients from quantum field theory |^6|, p7f . 

A different approach to the calculation of transport coefficients is the solu- 
tion of a Boltzmann equation (e.g. in relaxation time approximation) [^3| , |28f . 
However, this approach has the drawback of a zero width quasi-particle approx- 
imation at a very early stage of the computation, i.e., it is semi-classical. In 
contrast to such kinetic models, the formulation of the present work contains 
off-shell effects, and therefore the full quantum description, in a consistent way. 

Furthermore, some of the semi-classical calculations retain the relaxation 
time as an undetermined factor in the final expressions. While TFD in principle 
also provides an expression for the relaxation time (|, |(| , the method presented 
here is more direct and can be extended systematically. 

The new method was tested by calculating thermal conductivity and dif- 
fusion coefficient for the pionic component in a hot A-nucleon gas. Results 
were obtained, which are substantially higher than those found with the semi- 
classical zero width kinetic picture. As was pointed out in the previous section, 
part of this might be due to the too simplistic picture of the A-hole model with 
constant A-width. 

However, this finding is also consistent with an observation made, if the 
Kadanoff-Baym equations are solved directly rather than approximated by a 
Boltzmann equation: The relaxation of perturbations is slower, if the full quan- 
tum description is taken into account P9| . 
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Figure 1: Pion spectral function of the A- hole model at pb — 1.69 nuclear density 
and T = 0.1 GeV. 

Horizontal axis \k\ in GeV, vertical axis E in GeV. Plotted are lines of equal p(E, k) 
on 

logarithmic scale, see ref. [Jl8| for details. 
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Figure 2: Pion thermal conductivity x temperature. 
Thin straight lines: slopes oc T 3 , T 4 

Full thick line: A-hole model at pt, — 1.69 nuclear density; 
dashed line 0.92 and dash-dotted line 0.47 nuclear density. 

Dash-double-dotted line ref. Q, dotted line ref. Q, 
both from ir — it scattering data at T = 02g>i> = 0. 



Figure 3: Pion diffusion coefficient. 
Thin straight lines: slopes « T 3 , T 4 

Full thick line: A-hole model at p b = 1.69 nuclear density; 
dashed line 0.92 and dash-dotted line 0.47 nuclear density. 
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